Dynamic transition in driven vortices across the peak effect in superconductors 



O 
o 

(N 



Mahesh ChandranQ R. T. Scalettar, and G. T. Zimanyi 
Department of Physics, University of California, Davis, California 95616 
(Dated: February 2, 2008) 

We study the zero-temperature dynamic transition from the disordered flow to an ordered flow 
state in driven vortices in type-II superconductors. The transition current I p is marked by a sharp 
kink in the V(I) characteristic with a concomitant large increase in the defect concentration. On 
increasing magnetic field B, the I P (B) follows the behaviour of the critical current I C (B). Specifically, 
in the peak effect regime I P (B) increases rapidly along with I c . We also discuss the effect of varying 
disorder strength on I p . 
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PACS numbers: 74.25. Qt, 74.25. Sv 

The model of driven interacting particles in presence 
of quenched disorder captures essential features of the 
dynamics occurring in many condensed matter system. 
Of specific interest is the dynamics of vortices above 
the critical (or depinning) current I c in type II super- 
conductors. Early numerical simulation of the trans- 
port characteristics in 2D showed that the topological 
defects are generated close to I c , and are annealed at 
high driving force This suggested a dynamic tran- 
sition from the disordered flow to an ordered flow at a 
current I p > I c . Theoretically, Koshelev and Vinokur 
(KV) considered this transition at I p similar to an equi- 
librium "melting" transition of a clean system in the 
moving reference frame 0. The effect of quenched dis- 
order on the moving vortices is included in the theory 
through an effective shaking temperature T s h oc J . In 
Koshelev- Vinokur theory, the re-crystallization current 
I p increases as (T m — T) -1 on approaching the equilib- 
rium meltingtemperature T m of the static system. Later 
calculations 0,0,13 an d simulations d,0S0 identified 
the free flowing state above I p as the smectic phase. In 
the smectic phase, the particles move in static channels 
with quasi-long range order perpendicular to the flow 
compared with liquid like short range order within the 
channels |lfj. 

The dynamic transition in driven vortices was first 
studied experimentally in the transport measurements 
by Bhattacharya and Higgins (BH) [12j, and was later 
observed in neutron scattering [llj and Bitter decora- 
tion 0| experiments. BH studied the behaviour of tran- 
sition current I p as a function of the magnetic field B, 
particularly in the peak effect (PE) regime. The PE is 
marked by a rapid increase in I C (B), generally close to 
the upper critical field B&. BH found that the I P (B) in- 
creases rapidly in the field range in which the PE occurs. 
This behaviour of I P (B) in the PE regime is similar to 
the behaviour of I p (T) around T m as observed in simu- 
lation by KV in Ref. Q{. The similarity can be traced 
to the softening of the shear modulus cqq of the vortex 
lattice. But an important difference separates the two 
: unlike BH, KV do not see any corresponding increase 
in I C (T), i.e., no peak effect in the I C (T) is observed 
in the simulation. This makes it difficult to correlate 
the increase in I p (T) close to T m to the enhanced cou- 



pling of the vortex lattice to the quenched disorder. The 
increase in I P (T) could very well be due to increased 
thermal fluctuations on approaching T m , thus requiring 
larger currents to anneal the lattice defects. 

Recently, we showed that the PE occurs in a system 
of 2D vortices close to B C 2 at zero temperature [l4|. In 
this paper, we revisit the dynamic transition, particu- 
larly across the PE. We find that I p (B) increases rapidly 
in the field range in which the I C (B) shows the PE. The 
dynamic transition at I p is characterized by hysteresis in 
V(I) and a sharp peak in the dynamic resistance Rd{I). 
The topological defect concentration also shows a jump 
at Ip. We present the dynamic phase I-B diagram, and 
discuss the effect of increasing the disorder strength A 
on I p . 

Consider a 2D cross-section of a bulk type-II super- 
conductor perpendicular to the magnetic field B = Bz. 
Within London's approximation, the dynamics of vor- 
tices are governed by the overdamped equation of mo- 
tion 



dt 
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Here r) is the flux-flow viscosity, and the first term repre- 
sents the inter-vortex interaction given by the potential 

i 2 

U v {r) = 8 ^° X 'i K a (r/X), where K is the zeroth-order 
Bessel function, and f = (r 2 + 2^ 2 ) 1 / 2 . A and £ are the 
penetration depth and coherence length of the super- 
conductor, respectively. (f>o represents the flux quantum. 
The vortex pinning is added through the second term 
which is an attractive interaction with parabolic poten- 
tial wells U p (r) = Uo(~s — 1) f° r r < r p, and otherwise, 
centered at the random Rfc locations. The third term 
Fext — -J x <fio% is the Lorentz force due to transport 
current density J. The length is in units of A(B = 0) = Ao, 

,2 

and J is in units of c/o/0o where fo = gja^a ■ The time t 
is in units of 77A//0 whereas the velocity v of vortices is 
in units of fo jr\. 

We use the reduced magnetic field b = B/B C 2 with 
B C 2 = 2^p- 1 and is calculated from the lattice constant 
f = (4*) 3 (^)i The B dependence of A is follows 
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& 2 )2, with a similar expression for £. 
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The simulation is for the Ginzburg-Landau parameter 

A 



k = j = 10 which is typical of the low-T c materials 



The number of vortices N v was chosen in the range of 
800-1200. The prefactor Uq of the pinning potential is 
distributed randomly between A±0.01. The results pre- 
sented below are for the pin density n p — 2.315. We 
calculate the transport characteristics V(I), where the 
current I oc J y , and the voltage V oc (v x ). The topolog- 
ical defect concentration fd(I) is obtained by Delaunay 
triangulation of the real space configuration. We choose 
the free flowing vortex lattice at high driving currents 
as the initial configuration, and decrease the current to 
in small steps across I c . This particular method of 
preparing the system minimizes the influence of non- 
equilibrium defects which are present in field cooled sim- 
ulations. The equation of motion is time integrated by 
standard techniques. The V(I) and fd(I) are time aver- 
aged in the steady state. Parallel algorithms were em- 
ployed to speed up the simulation at high densities . 



0.4 




1 










0.2 






- 




W...1......... 



(a) 




FIG. 1: (a) The V(I) characteristics for A = 0.04 and 6=0.7. 
Inset : The defect concentration fd(I), and the dynamic re- 
sistance Rd{I). (b) Hysteresis in V(I) and fd(I) (upper in- 
set) across I p . The lower inset shows the {8y c h n i) / a-a , where 
Qq = ^-Oo and ciq is the lattice constant. 



between I C <I <I p where some vortices remain immobile 
and large transverse excursions of the active channels are 
present, and (c) an ordered flow state for I>I p with all 
vortices moving and interchannel hopping takes place 
with, at most, neighboring channels. The fd(I), and the 
dynamic resistance Rd(I) = is plotted in the inset of 
Fig. 1(a). The current I p is marked by a sharp peak in 
Rd{I) and a large change in fd- The response function 
V(I) shows a kink at I p , as evident from Fig. 1(a). We 
find that V(T) and fd(I) are also hysteretic across I p , 
as shown in Fig. 1(b). The kink and the hysteresis in 
the response function V(I) suggests that the dynamic 
transition occurs at I p . The transition is of first order 
nature, as was identified in Ref. 0. 

In Fig. 2, real space configurations are shown for 
currents above and below I p corresponding to param- 
eters in Fig. 1. The instantaneous configuration shows 
a relatively ordered lattice above I p . The vortex tra- 
jectories are well separated channels aligned parallel to 
v x . The average transverse wandering of the channel 
(Sy c hni) < olOq, where a$ — ^a-o and ciq is the lattice 
constant, and a~0.05. Just below I p , a becomes « 0.3, 
as shown in the inset of Fig. 1(b). The slowing down of 
vortices effectively couples the transverse velocity com- 
ponent v y to the quenched disorder. This induces large 
scale dislocations, and the dynamic friction (oc Rd) in- 
creases. The increase in transverse wandering of vortices 
is also reflected in the dense braiding of the active chan- 
nels. The behaviour of (Sy c hni) across I p is similar to 
the Lindemann criterion for thermal melting. Note that 
some of the vortices appear as immobile for / = 0.055 
in Fig. 2(b) (rightside panel). The region of immobile 
vortices grow with decreasing / until I c is reached when 
the whole system is pinned. 




FIG. 2: Instantaneous vortex configuration (left) and vortex 



trajectories (right) for (a) I- 
The 6 = 0.7, and A = 0.04. 



0.06>/ p , and (b) 7 = 0.055<7 P . 



Fig. 1(a) shows the V(I) curve for b = 0.7 and A = 
0.04. Broadly, three current regimes can be identified: 
(a) a pinned state for I<I C , (b) a disordered flow state 



The Rd{I) has been us ed p reviously to identify the 
phase boundary. In Ref. [12J, the transition from the 
disordered flow to the ordered flow state is identified as 
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the current at which Rd(I) ~ Rbs where Rbs is the 
asymptotic Bardeen-Stephen flux-flow resistance. This 
current is marked as I cr in Fig. 1(a). In Ref. 8], the state 
between I p and I cr is identified as the smectic phase, 
whereas above I cr , the state is defined as a transversely 
frozen phase. The later state is distinguished from the 
smectic phase by the absence of transverse jumps by vor- 
tices between adjacent channels. 

We note that neither V(I) nor fd(I) shows any feature 
at I cr . Therefore, a more appropriate description of the 
state between I p and I cr is possible in terms of fluctua- 
tion which allows occasional transverse excursions of the 
vortices between adjacent channels. The rate of such 
transverse jumps can be calculated by employing Ar- 
rhenius relation which is conventionally used to describe 
thermal activation over energy barrier, but with temper- 
ature replaced by T s h [l^. Such a description can ac- 
count for the suppression of transverse jumps above 7 cr , 
and also for the experimentally observed Rd(I) curves. 
To summarize, the dynamic transition in the moving lat- 
tice occurs at I p at which the dynamic resistance Rd(I) 
shows a peak. The flowing state above I p is a smectic 
phase with small longitudinal correlation and quasi-long 
range order in the transverse direction. 
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FIG. 3: V(I) curves for A = 0.04 for fields across the PE. 
Inset : V(7) curves for increasing disorder strength A for 
b = 0.6. The curves in the inset are shifted horizontally by 
0.01. The arrows mark the transition current I p . 

We now consider the behaviour of I p as the magnetic 
field b = B/B C 2 is changed, particularly across the PE in 
I c (b). For A = 0.04, the onset field for PE is b op ^0.75 
whereas the peak occurs at b pea k ~ 0.9. Fig. 3 shows 
the V(I) curves for four values of b between 0.5 and 0.8. 
The I p can be easily identified by the sharp kink in the 
V(7) curves. With increasing b, the I p first decreases 
before increasing rapidly in the PE field regime. The 
fd(I) and Rd(I) characteristics across I p are similar to 
that shown in Fig. 1 for all values of b. In the inset of 
Fig. 3, the effect of increasing disorder strength A on I p 
is shown for b = 0.6. Other than the expected increase 
in I p with increasing A, we find that the transition at I p 



is broadened for A > 0.05. Also, the hysteresis in V(I) 
across I p is too small to be distinguished. This suggests 
that the first-order nature of the dynamic transition at 
I p does not survive for A > 0.05. 

Fig. 4(a) shows the dynamic phase diagram I-B for 
A = 0.04 which summarizes the main result of the paper. 
The I p (b) and I c {b) shows similar behaviour as b ap- 
proaches the upper critical field value. Particularly, the 
rapid increase in I p (b) coincides with the peak effect in 
7 C (6). The region between I p (b) and I c (b) constitutes the 
disordered flow regime, whereas above I p vortices flow in 
ordered channels. For A = 0.04, the I p (b) = fc(6 - b)^ 1 
with bo = 0.92 and k = 0.015 gives a reasonable fit for 
b>b op , as shown by the thick line in Fig. 4(a). This form 
of I p (b) was motivated by the behaviour of I P (T) on ap- 
proaching the equilibrium melting transition in Ref. |2|. 
The I-B plot in Fig. 4(a) thus supports the picture that 
both the static and the dynamic friction of the vortex sys- 
tem increases rapidly in the field range in which the peak 
effect occurs. As shown in Ref. 01 > the PE in the critical 
current I c is driven by the softening of the vortex inter- 
action on approaching the upper critical field B C 2- Soft- 
ening of the vortex interaction also reduces the shaking 
temperature T s h ocl^ 1 required for the dynamic ordering 
(or recrystallization) of the vortices which explains the 
increase in I p (b) in the PE regime. Thus, an increase in 
I c (as a function of B or T) implies an increase in I p . We 
emphasize that the reverse does not hold, i.e., increase 
in Ip does not imply an increase in I c . This can be seen 
from the behaviour of I P (T) and 7 C (T) in Ref. Q where 
thermal depinning causes I C (T) to decrease monotoni- 
cally even as I P (T) increases on approaching the equi- 
librium melting temperature T m . Overall, we find the 
behaviour of I P {B) and I C (B) in good agreement with 
the experimental observation [l2j. 

Fig. 4(b) shows the I c and the I p as a function of 
A for b = 0.6. For the system size used in the simula- 
tion (N v ^1000), the depinning transition is elastic for 
A < 0.03 (hence, I p = 7 C ), whereas for A > 0.05, I p cx A. 
The two disorder regimes can be identified as the weak 
pinning and the strong pinning regime, respectively. The 
crossover from weak to strong pinning occurs for A be- 
tween 0.03 and 0.05, and is also the range in which a 
sharp dynamic transition is observed. We find that for 
the intermediate pinning the real space configuration at 
1 = show domains of ordered lattice separated by do- 
main walls. The presence of ordered region ensures that 
the distribution of v x oc is narrow which consequently 
gives a sharp dynamic transition. On the other hand, 
increasing A above 0.05 decreases the domain size to ~ 
2-3ao. This leads to a broader distribution of v x which 
in turn broadens the dynamic transition. Thus, a sharp 
dynamic transition at I p implies the existence of ordered 
regions with a narrow distribution of the domain size. In 
Ref. ^3 > the coexistence of the ordered and the disor- 
dered regions was shown to underlie the shape of the 
V(7) curves in the vicinity of the PE. The disordered 



4 




(b) B/B c2 = 0.60 




0.12 



FIG. 4: (a)The dynamic phase diagram showing I p (b) and 
I c {b) for A = 0.04. The thick line is a fit to I p (b) ~ (6 - b)' 1 
in the peak effect regime, (b) The I p and I c as a function of 
A for 6 = 0.60. The dashed line is a linear fit to the data in 
the range A > 0.05. 

state in Ref. 01 is thought to appear due to the injec- 
tion of vortices across the surface barrier. Our simula- 
tion suggests that the shape of the V(I) curve depends 
crucially on the distribution of the size of the ordered 
domains in equilibrium. 

In conclusion, we have shown that for disordered type- 
II superconductors, the dynamic transition current I p 
and the critical current I c shows similar behaviour as the 
magnetic field B is varied. I p decreases with B for fields 
below the PE regime. In the peak effect regime, I P (B) 
increases rapidly along with the I C (B). The dynamic 



transition is sharp for the intermediate pinning strength 
due to the presence of large domains of ordered vortex 
lattice in equilibrium. 
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